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I. INTRODUCTION 

The nature of dark matter and dark energy is one of the most important issues today in physics. 
There are strong observational evidences in astrophysics and cosmology for the existence of these 
two exotic components of the cosmic energy budget, indicating that about 95% of the universe 
is composed of dark matter (about 25%) and dark energy (about 70%), but no direct detection 
has been reported until now. The usual candidates to dark matter (neutralinos and axions, for 
example) and dark energy (cosmological constant, quintessence, etc.) lead to very robust scenarios, 
but at same time they must face theoretical and observational issues. For recent reviews on the 



subject, see for example 



3] 



Among the alternatives to the standard description of the dark sector there is the possibility of a 
modification of gravity theory on large scales, see J] for a recent review on the subject. An example 
are the so-called f (R) theories, which are based on the inclusion of non-linear curvature terms in 

n 

the Einstein-Hilbert action [5jJ. Another example, the one we pursue in this paper, is to touch 
one of the cornerstone of the gravity theory: the usual conservation laws for matter components. 
This kind of formulation was introduced by Rastall some 40 years ago [f], Q], and has been recently 
investigated in a cosmological context, giving some interesting results concerning the dynamics of 
the dark sector [s-12]. 



RastalFs motivation for modifying the usual conservation laws is based on the fact that the 
latter have been directly tested only locally and in a weak-field regime. On the other hand, the 
introduction of covariant derivatives imply, in some sense, an exchange of energy between matter 
and the gravitational field. Hence, in general, non-trivial generalizations of the conservation law 
are in principle possible. Besides, particle production in curved space-times is a central issue 
of quantum field theory on such spaces. Therefore, one may regard modifications of the usual 
(classical) conservation laws as effective, semi-classical approaches to such phenomenon. RastalFs 
proposal is the following: 

= kR> v , (1) 

where the semicolon denotes the covariant derivative with respect to a generic metric and k is a 
(dimension-full) free parameter. The above relation can be rewritten as 

= ^1"'" > (2) 

where T is the trace of the stress-energy tensor and 7 is now a dimensionless free parameter. When 
K = 0, then 7 = 1 and the usual conservation law (and thus General Relativity) is recovered. 
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In a one-fluid model, it is possible to redefine the energy-momentum tensor in order to recover 



the usual conservation law 13j]. In this sense, Rast all's theory is just a redefinition of the fluid 
equation of state. However, in a multi-fluid case the modification introduced by Rastall opens 
possibilities of non-trivial interaction between the different components. In [9J], this quality is used 
to investigate the consequences of the modified conservation law for a model of the dark sector of 
the universe. The authors investigate a two-fluid model, one of them being pressure-less matter 
Pm = 0, whereas the other obeying the vacuum energy equation of state p x = —p x . Assuming that 
the matter component obeys the usual conservation law, then the vacuum energy conservation law 
is affected by the presence of matter, via Eq. ([2]). The main result of [9] indicate that the model 
is completely equivalent to the ACDM at the background and linear perturbations levels. There 
is just one striking difference: dark energy may now agglomerate. This fact could have important 
consequences at the non-linear level, which is a regime where the ACDM faces some difficulties 

Q. 

In the model studied in the equivalence with the ACDM at background and linear per- 
turbations levels implies that no constraints on the parameter 7 can be established using the 
corresponding observational tests. Such constraints are, on the other hand, possible using non- 
. inear data. In the present paper our goal is to verify to what extent the model studied in reference 
9] is a favourable configuration. To do this, we repeat the analysis made there but with the x 
component now being described by a more general equation state, i.e. p x = oo x p x , with lo x constant. 

We show that in this case the parameter 7 appears explicitly in the background and linear per- 
turbation equations. Therefore, using type la supernovae and matter power spectrum observational 
data, we calculate its probability distribution function (PDF) along with the one for the matter 
density parameter O m o> and for the equation of state parameter u x . We obtain that the PDF's 
are are peaked about u) x ~ —1, £l m o ~ 0.3 and 7 ~ 2. Remarkably, the case 7 = 2 corresponds to 
the original formulation of General Relativity (based on the Ricci tensor) and also to a zero speed 
of sound for a self-interacting scalar field obeying the modified conservation equation proposed by 
Rastall's theory lo| . 



The paper is organized as follows. In Sec. IH we present Rastall's theory, deriving some cosmo- 
logical considerations. In Sec. lIIII the fields equation are presented and in Sec. HVI the confrontation 
with observation is carried out. In Sec. [V] we present our conclusions. 



4 



II. RASTALL'S THEORY 



According to Eq. ([2]), Einstein equations must be modified as follows: 



Guv — Rfiv — -^QnuR — 8ttG ^Tfj,v — — dfiuT^j , 



(3) 



V 2 

in order to be compatible with Bianchi identities (we use c = 1 units hereafter). 

We consider a two-fluid model. The first component mimics baryons and dark matter, i.e. has 
negligible effective pressure, whereas the second describes an exotic dark component responsible 
for the acceleration of the universe, and has an equation of state of the form p x = co x p x , with uj x 
constant. We assume the matter component to conserve as usual. This is important in order to 
have matter agglomeration required to form the local structures. The field equations become 

1 



R»v = 8vrG 



^ + r---(2- 7 )^(r- + r m ) 



1 x ;m 



7-1 



V^x ' J -mJ i ± m ;M — u ' 



(4) 
(5) 



where subscripts or superscripts x and m denote the dark energy and the matter component, 
respectively. On large scales we assume the universe to be described by the Friedmann-Lemaitre- 
Robertson- Walker (FLRW) metric, 

dr 2 



ds z 



dt 2 - a(ty 



+ r z (d9 z + sin z 9dip z 



(6) 



1 — kr 2 

where a(t) is the scale factor and k is the curvature of the spatial sections. The perfect fluid 
energy-momentum tensor has the following form: 



(7) 



where j = x, m. We focus on the flat case (k = 0), where the field equations take on the following 
form: 

a 4?rG 



3 

4ttG 



{[ 3 (7 - 2 )W X ~ 7] Px ~ IPm} , 

Px [3 - 7 - 3(1 - j)w x + (3 - j)p r , 



d 'y — 1 

p x + 3-(l + w x )p x = —— {(1 - 3w x )p x + p m } 
a 2 



i a 

Pm ~T '~>~Pn 

a 







(8) 
(9) 
(10) 
(11) 



where the dot denotes derivative with respect to the cosmic time t. Note how the expansion history, 
therefore the accelerated phase, of the universe depends of the parameters 7 and oj x . The evolution 
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profiles of p m and p x are easily determined by solving the last two equations, yielding 

6(1 + ^) 3d— ~/)0rr,n Q 

p x = Px0a -*-«-w-*°.> + - /J- 1 7^7° - . a~ 3 , (12) 

6w x + 3(7 - 1)(1 - 3w x ) 

PmO /.. q n 

= — a - , (13) 

a 6 

where p m o and p^o are their present values (i.e. for a = 1). The present-time density parameters 

of the exotic fluid x and of matter are given by 

8irGp x0 _ 8irGp m 

q ' m0 ~ 3# 2 

and in a flat universe are related as follows 



&x0 - , ^mo - 5 ( 14 ) 



2-0 



6^+3(7 V(i-3^) [3 - 7 " 3(1 - 7R] + (3 - 7)} 



= ^ ^r: ; • (15) 

3 - 7 - 3(1 - 7)w x 

In the special case where 7 = 1 the above relation recovers the corresponding General Relativity 
one f2 x o + ^mo = 1- We now turn to the perturbation analysis of our two- fluid model. 

III. PERTURB ATI VE ANALYSIS 

We write the perturbed metric as = gffl + h^ u , where gffi indicates the background flat 
FLRW metric and hav is a small fluctuation. We choose the synchronous gauge condition, i.e. 
hfj,o = and introduce the perturbations in the matter functions as follows: 

Px = Pi 0) + ?>Px , Pm = Pin! + SPm , U m = U$ + 5u m , (16) 
Px = Pi 0) + $Px , Pm = P$ + SPm , Ux = m! 0) + $U X , (17) 

In the expressions (fT6|) - (fT7]) . the superscript (0) indicates the background functions; 5p x , 5p m , 5u m , 
5u x , 5p x , 5p m , represent the perturbed quantities in density, four- velocity and pressure. We also 
introduce the following customary definitions: 

S x = , 6 m = , 9 = d k 5u k , h = — . (18) 

Px Pm a 

For deeper detail, see for example [15]. After standard calculations, the perturbed conservation 
equation for the matter component can be cast in the following very simple form: 

L = \ ■ (19) 

Introducing the functions 

5 (a) = a = -aj[ 7 + 3(2-7K]^ + 7^4 , (20) 



/( „ ) a Jtk[ 3 - 7 - 3(1 - 7)0;,] + ^(3 - 7) , (21) 
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where 



Q x = Q x0 a 2-(7-l)(l-3^) + 



3(1-7)0 



mO 



6w x + 3(7 - 1)(1 - 3w x ) 



(22) 
(23) 



and calculating the perturbed component i?oo> we get the following equation: 

J (ft) ■ ^'T 3 
5 m + 2 5 m — — O m <5 m = -[7 + 3(2 - 7)w ;c ](5 :!; :c . 

Calculating the perturbation of ([5]) we obtain for v = 



(24) 



(5^ + 3^-^(1 + u) x )6p x + (1 + a; a .) / o a: ( G x - ~ 
a V 2 



7-1 



and for ^ = i 



(1 + u x )p x 9 x + (1 + u^)/^©* + (1 + w £C )p a; ( 1 (-), 



[(1 - 3uj x )5p x + 6p m ] 



•/(«) 



(25) 



V 2 ^ 



1 — 7 / 37 — 5 

+ I — ~ — I w x 



+ 



V 2 5p m ( 1 - 7 



(26) 



It is useful to rewrite the above equations in terms of the derivative of the scale factor a instead 
of the comic time t since the former is directly connected with the dimensionless redshift quantity 
through z = — 1 + 1/a. In this way, Eq. (I24h becomes 



ft' + ft ( g ( Q ) 1 2 



' m " m V/ 2 («) a/ 2/2 (a) " m 2/2(0) 
Equation (|25|) turns into 



[7 + 3(2 - j)u x ] , 



(27) 



51 



1 



1 + ^(1-3^ 



-3^(1 + 



<5L 



(l + ^) + lfe(7 



1 + 



(l-7)(l-3^) 



/(« 



+ 



+ 



1+^2 



<M7 - l)^r 



2/(a)^ 



1 + 



(l-7)(l-3c^) 



(28) 



and finally Eq. (126ft reads 



©' + 



Px a 



e. 



1 v 2 ^ 

1 + u x f(a)a 2 



1 — 7 _l_ ^37 — 5 



w 3 



+ 



1 O m V 2 <5 m / 1 - 7 



1 + ft x f(a)a 2 \ 2 



(29) 
(30) 



where the prime denotes derivation with respect to the scale factor. 
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IV. OBSERVATIONAL CONSTRAINTS 

We submit our model to two tests: the first focuses on the background evolution and is based 
on type la supernova data (the Constitution set [la]): the second is a perturbation test and we 



employ matter power spectrum data from the 2dFGRS 



17] 



We perform a Bayesian analysis, calculating first the x 2 function, defined as follows: 

x 2 = E { 1 J l) ■ ( 31 ) 
i=i 1 

where ./V is the total number of observational data, xf 1 is the theoretical prediction of some physical 
quantity and xf is its value observed within an error bar erj. Assuming the data to be independent 
Gaussian random variables, see for example [3], the likelihood distribution function is constructed 
from the \ 2 function as follows 

L(p>) = e -xV)/2 , (32) 

where p 3 represent the free parameters of the theory. For our model, we have four free parameters, 
i.e. h, £l m o, 7 and u x . However, we choose to fix h = 0.7, on the basis of the current consensus on 



the value of the Hubble constant, see also 
space to three. 



191 ] , hereby reducing the dimension of the free-parameter 



A. Supernovae la 

Type la supernovae data consist in the distance modulus fi, i.e. 

^ = m obs ^i)-M = 5log(^\ + 25, (33) 

where m i, s is the apparent magnitude, M is the absolute magnitude and cLl is the luminosity 
distance. It can also be written as 

H = 5log 10 D L (z) + no , (34) 
where Dl = (Hq(1l)/c is the Hubble-free luminosity distance and is the zero point offset, defined 

by 

Mo = 5 log 10 (j^A + 25 = 42.38 - 5 log 10 h . (35) 

We employ data from the so-called Constitution set [ie|, which includes 397 distance moduli, out 
of which 100 come from the new low-z CfA3 sample and the rest from the Union set. Both samples 
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have a redshift range of 0.015 < z < 1.55. The main improvement of the Constitution sample 
is the inclusion of a larger number of nearby (z < 0.2) supernovae, thus reducing the statistical 
uncertainty. 

The chi-squared we employ for the type la supernovae test is then 

397 2 

*Wr) = E '"*■"'' *> ^ , (36) 

— C ■ • 

1=1 ofe,i 

where p = (ri m o, Wj, 7, /io)- The chi-squared can be minimized with respect to fj,o, since the latter 
is an independent parameter. Expanding Eq. (f36|) with respect to //q, we obtain 



X 2 (pW = A(p) - 2 Mo B(p) + ^Cfr) , (37) 
which has a minimum for ^0 = B(p)/C(p), giving thus 

XSNIa,rnin = XSNIa = ^(p) ' 

where 

A(pi) = y i^(p^ = o)-^] 2 (39) 

= y = q)-»<*. (40) 

C{pi) = \ . (41) 

i 

As a side remark, for the single-fluid model Rastall's theory has been confronted against supernova 
data (Union sample) leading to results that can be competitive to the ACDM data under special 
conditions [^(J. 



B. Power Spectrum 

We use the 2dFGRS data [17] to compare the predictions on the matter power spectrum of our 
model with observation and obtain constraints on the parameters. The chi-squared, in this case, is 

Xp( k )^ l m0,^x,'l) = 2 , 2 • \ A2 ) 

t=l a i 

The number of free parameters is three (co x , 7 and Q m o) because h is absorbed in the units of the 
wave-number k (h Mpc -1 ). 
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In order to compute the matter power spectrum, we choose the BBKS fitting function as initial 
conditio, after ecuant, ft, the tracer function ft 

rp (k s = ln(l + 2.34g) 

bbks I ) 2 .34g[l + 3.89 9 + (16.1 g ) 2 + (5.46 9 ) 3 + (6.71g) 4 ] 1 /4 ' 1 > 

where q = ^/(E/iMpc -1 ), with the shape parameter £ = Sl m o^ e ^ M properly taking into account 
the baryon contribution, which we assume Sl^o = 0.04. The matter power spectrum can be finally 
written in the form 

P(k) = Ak n T 2 (k)T 2 BKS (k), (44) 



19, 



23|), and 



where n measures the slope of the primordial power spectrum (we assume n = 1 
A is a normalization constant. To obtain power spectra for our model, we numerically solve the 
relevant equations presented in the previous section from z = 500 to z ^ and calculate T 2 of 
Eq. (|44p . For the normalization constant, we follow the methodology of 



Note that in 



23] 



the matter power spectrum is used to constrain the parameter 7, but with a 



decomposition of the matter components different from that employed here and in 
in 7 = 1, i.e. General Relativity, as the preferred configuration. 



This results 



C. Joint analysis 

Considering the observational data to be independent, we define the total chi-squared as 

X 2 = XsNIa + Xp( k ) > ( 45 ) 

and the likelihood as 

L(n m0 , u x ,j) = exp (-x 2 /2) • (46) 

We assume flat priors for the parameters, therefore, up to a normalization constant, the likelihood 
is also the posterior probability. In FIG. [THZ1 we show the results for the individual and joint 
analysis. The principal result is the low probability for the fi m o parameter using the 2dFGRS data 
only and a peak of probability near 7 ~ 2 and lo x ~ —1. 
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Data 


Xmin ^m0 






7 


2dFGRS (7 = 2.0) 


18.1400 0.200^;5^ 


-0.955^ 


-0.120 

-0.031 


2.0 


SNIa (7 = 2.0) 


465.508 0.310±g:|25 


-0.985^ 


-0.150 
-0.120 


2.0 


2dFGRS+SNIa (7 = 2.0) 483.657 0.30ltg;^g 


-l.OOSj 


-0.130 
-0.120 


2.0 


2dFGRS 


18.1497 0.239±g;£j$ 


-1.350j 


-0.107 
-0.205 


1.810l Jo7 


SNIa 


465.294 0.306±g Jig 


-0.965^ 


-0.201 
-0.051 


1 001 +0.305 


2dFGRS+SNIa 


482.649 0-21±g;g|| 




1 01 4+1.205 
i - 8i4 -1.502 



TABLE I. Observational constraints on the model parameters (la errors). 



V. RESULTS AND DISCUSSION 



We compared the predictions of a cosmological model based on Rastall's theory against type la 
supernovae and matter power spectrum data. We determined a probability distribution function 
for each of the free parameters (£l m o, u x , 7). 

In Table 1 we show the best fit values (with la errors). They are very close to the ACDM 
results, but with a quite important dispersion. A clear feature is that supernovae test favours a 
higher value for f2 m o with respect to the power spectrum estimations. The joint analysis, for this 
parameter, depends on fixing or not 7: it is higher when 7 = 2 and smaller when 7 is left free. 
It is clear also from table 1 that the PDF for 7 is not highly dependent on which test is used. 
In general, 7 is peaked around 1.8 but with a quite large dispersion, which leads to the value 
7 = 2 quite probable. Remarkably, the General Relativity case 7 = 1 is not statistically favoured 
from the analysis performed here. For the equation of state parameter in almost all cases, the 
most probable configuration is near the cosmological constant case oj x = —1, except when only the 
power spectrum test is used, with 7 free: for this case, the equation of state parameter extends 
considerably to the phantom region (u x < — 1). 

The results reported here indicate a scenario that fits the requirements for the structure forma- 
tion process, since 7 = 2 implies a sound speed around zero lOJ], at least for the scalar field case. 
Moreover, the vacuum energy equation of state is somehow favoured, supporting the configuration 
studied in [9]. In this sense, a natural extension is to consider the non-linear regime for structure 
formation. Moreover, Rastall's theory has been investigated until now on a cosmological ground 
only. It is important to verify its properties on smaller scales, for example investigating the stability 
properties of objects such as stars and the weak-field dynamics. 
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FIG. 2. Marginalized probability density functions for the parameters u x and fi m o for fixed 7 = 2 and using 
type la supernovae data only. 




n m „ y 



FIG. 3. Marginalized probability density functions for the parameters uj x and £l m o and 7, using type la 
supernovae data only. 
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FIG. 4. Marginalized probability density functions for the parameters uj x and fi m o for fixed 7 = 2 and using 
matter power spectrum data only. 




y a„„ 



FIG. 5. Marginalized probability density functions for the parameters uj x and f2 m o and 7, using matter 
power spectrum data only. 




FIG. 6. Confidence contours at la and 2a levels on the (f2 m o> uj x o) plane and marginalized probability density 
functions for the parameters w x and fi m o- I n these plots we fix 7 = 2 and use both type la supernovae and 
matter power spectrum data. 




FIG. 7. Confidence contours at la and 2a levels on the (£l m o,w x o) plane and marginalized probability 
density functions for the parameters u x , ft m o and 7 = 2. We use both type la supernovae and matter power 
spectrum data. 
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